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Abstract 

We present two integrable discretisations of a general differential-difference 
bicomponent Volterra system. The results are obtained by discretising di¬ 
rectly the corresponding Hirota bilinear equations in two different ways. 
Multisoliton solutions are presented together with a new discrete form of 
Lotka-Volterra equation obtained by an alternative bilinearisation. 

1 Introduction 

One of the main difficulty in the topic of integrable systems is obtaining an 
integrable discretisation of a given partial differential or differential-difference 
integrable system. Applying integrability criteria like complexity growth [Tj or 
singularity confinement [2] is not always quite easy because the discrete lattice 
equations have, in general, complicated forms. However, is much more con¬ 
venient to start with some general lattice equation and impose some simpler 
(but more restrictive) integrability requirements (like, for instance, cube con¬ 
sistency [3] which leads immediately to the discrete variant of “zero curvature 
representation”). 

Among the methods of finding integrable discretisations, one of the powerful 
one is the Hirota bilinear method. The idea is quite simple. First, the integrable 
differential or differential-difference integrable system has to be correctly bilin- 
earised (in the sense of allowing construction of general multisoliton solution). 
After that, in a first step, one has to replace differential Hirota operators with 
discrete ones preserving gauge invariance. Of course the resulting bilinear fully 
discrete system is not necesarily integrable, so in the second step, the multisoli¬ 
ton solution must be found [3] . If this exists then the discrete bilinear system is 
integrable and, in the final step which is rather complicated, the nonlinear form 
has to be recovered. Here, it is possible to introduce an auxiliary function, as 
Hirota have shown in [5j. We have to clarify here the concept of integrability. 
Usually, the integrability for a given partial differential or partial discrete equa¬ 
tion, is relying on the existence of an infinite number of independent integrals 
in involution that usually can be computed from the Lax pairs. Accordingly, if 
the equation can be written as a compatibility condition of two nontrivial linear 
operators (i.e. Lax pairs) then it is automatically completely integrable. How¬ 
ever, the alternative formulation given by Hirota bilinear form, is also used in 
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proving integrability by requiring the existence of general multisoliton solution. 
Here the word general means the solution describing multiple collisions of an 
arbitrary number of solitons having arbitrary parameters and phases and also for 
all branches of dispersion relations. Any constraint in the form of multisoliton 
solution breaks complete integrability. The integrals of motion are related to the 
soliton parameters which, because of elastic collisions, remain unchanged (and 
they are also related to the spectrum of Lax operator invariant as well since the 
equation itself is an isospectral deformation). Construction of the multisoliton 
solution is quite difficult, but it has been observed that once three general soli¬ 
ton solution is constructed then it can be proved by induction that the general 
N-soliton can be constructed as well. The existence of three-soliton solution has 
been used in classification of completely integrable bilinear equations m, m 
as an integrability criterion. 

In this paper we are going to give two integrable discretisations of a gen¬ 
eral two component bidirectional Volterra system (in the sense that the dis¬ 
persion relation has two branches). This system is rather old. It has been 
formulated for the first time by Hirota and Satsuma [7] but with constraints on 
parameters. Then, various solutions (rational, white and dark solitonic) have 
been obtained 0 , 0 , m- Curious enough, the fully discretisation has not 
been studied so far, although the initial system can be formulated as a cou¬ 
pled differential-difference focusing or defocusing mKdV system with non-zero 
boundary conditions. In this paper we present two fully discrete bilinear forms, 
show the structure of V-soliton solution and then recover the nonlinear form. 
For the second discretisation we extend the approach of Hirota presented in j5] 
and use two auxiliary functions. An important fact is that the multisoliton so¬ 
lution in the case of second discretisation has practically the same phase factors 
and interaction terms as in the differential-difference one. The same fact has 
been observed by Hirota and Tsujimoto in mm, ® for other cases. They 
have shown that for many examples including lattice mKdV, lattice NLS, lat¬ 
tice coupled mKdV, the structure of soliton solution remains the same as in the 
case of differential-difference analogs (although, it is true that all their examples 
support only unidirectional solitons). We think that this fact happens in many 
other cases and as an other example we construct a new completely integrable 
discretisation of Lotka-Volterra equation. 

2 Discretisation of general Volterra system 

The differential-difference system under consideration is: 

Qn = (^0 T Cl Q n T C2Q n ){Rn+l Rn— l) (1) 

Rn = (cq T ClRn T C2-^ n )(^n+l Qn— l) 
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where Co, Ci, C 2 are arbitrary constants and the dot means derivative with respect 
to time. In [8] we have shown that following the simple scalings and translations, 


2c 2 , cr _ 2 c 2 

U n — r- - -(Qn + )j V„ — /- - j 

v 4c 0 c 2 - cf 2c 2 yj 4c 0 c 2 - c{ 


(Rn 


Cl , 
2c 2 ' 


we can cast the system in the following form: 


— (1 H" 'U J n)(Vn -\-1 V n — l) (2) 

V n = (1 + Vn)(u n + 1 - «n-l) 

where u n ,v n -4 a as n —4 ±oo and a = ci(4coc 2 — c 2 ) -1 ' 2 . These transforma¬ 
tions are valid only in the case of nonzero c 2 . Of course, one must be careful 
about the square root appearing in the definition of a. In case of negative 
argument the system is changed in a defocusing form: 


tin — (I 1 'Cn— l) (3) 

tin — (1 'C n ){u n -\-1 Un— l) 


From the point of view of integrability this is not crucial since an overall imag¬ 
inary unit factor for u n , v n will change the system one into the other. But for 
the real soliton dynamics (which we are not discussing here)the solutions of the 
above systems are behaving completely different (only the unicomponent case 


has been analysed quite recently in [lU]h 

For the system (2) the Hirota bilinear form is given by: 

fnQn fn9n — (1 T Cl ) (Fn+1 G n —\ G n -\-lF n —i) (4) 

F n G n ~ F n G n = (1 + a 2 )(/„+i<?n-l — gn+lfn-l) (5) 

(1 + ia)F n+ iG n -i + (1 — ia)F n -iG n +i = 2f n g n (6) 

(1 T icy) f n +ig n —i T (1 iot) /n—i<M+i = 2 F n G n (7) 

where the nonlinear substitutions are: 


1 d f n {t) i d F n {t) 

tin = Cl — — — 111 7TT, v n = a - - — In 

2 at g n (t) 2 ot G n (t) 

The bilinear system (4)-(7) is completely integrable 0 and has the following 
IV-soliton solution: 
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where the dispersion relation and phase factors are: 
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e iPj V 
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In order to construct an integrable discretisation we replace time derivatives 
in (4) and (5) with finite differences (t —> m) 


L -t ^ ( f(n , m + d) - f(n, m)) 

and impose the invariance of the resulting bilinear equation with respect to 
multiplication with exp (pun F vm ) for any (bilinear gauge invariance). In 
this first discretisation we discretise also the equations (6) and (7) by assuming 
a gauge-invariant shift in the time variable. 

The fully discrete gauge invariant bilinear equations are given by: 


/n<?n — /n<?n — <5(1 + O 2 ) {F n +1 G n -1 — G n +lF n -i) (8) 

F n G n — F n G n = 5(1 + a 2 )(/ n+ ig n _i — g n +ifn-i) (9) 

(1 F iot)F n +\G n —\ T (1 ia)F n _ 1 G n+1 = fnQn F fn9n (19) 

(1 F lot') T (1 icx'yfn— F n G n -)- F n G n (H) 


where f n = f(n,m F S ) etc. In order to check the integrability one has to 
compute the general 7V-soliton solution of the above system. With the aid of 
a symbolic computation software like MATHEMATICA one can easily find 3- 
soliton solution. It can be extended to N-soliton solution in the form (the proof 
is given in the Appendix): 
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where the dispersion relation and phase factors are given by: 

Pi + (5ei(l + a 2 ) 


<h = 


Pi + Pi £ i^(l + « 2 


(14) 

(15) 


(16) 


ia{-l + ^( Pl + Pi 1 )) ej ia(-l + f (pi +Pj *)) £j 

-p) -1 )^ + 6 +(5a 2 ) 2 ’ * (Pi -Pj _1 )(l + <5 + <5a 2 ) 2 
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Accordingly, our bilinear system is an integrable one. Now we can proceed to 
recover the nonlinear form. Dividing (8) by (10) and (9) by (11) and taking into 
account that tan(| log (G/F)) = i(G — F)/(G + F) for any G and F we obtain 
the following system: 


tan(Q n Qn) 
tan(l? n Rn) 


<5(1 + a 2 ) tan(E ra+ i - R n -i) 
1 + atan(^ ra+ i - R n -i ) 

(5(1 + a 2 ) tan(Q w+ i - Q n _ i) 
1 + atan(Q n+1 - Qn- 1) 


(19) 


where Q n = | log(/„/p n ), R n = | log (F n /G n ). To our knowledge this system is 
a new one. We do not know how it is related to other integrable discretizations 
of Volterra systems. In the case of a —> 0 the classical lattice self-dual network 
of Hirota is obtained. 

Since the phase factors are defined up to the multiplication with the same 
constant factor (which we take to be the imaginary unit) then f n , g n (and F n . 
G n as well) will be complex conjugated so the physical fields Q n and R n are 
real functions. Here we give the 1 and 2-soliton solution (see also the figure): 
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Figure 1: One and two-soliton solution for Q n ( p\ = 2,p 2 = 0.7, a = 1 , 5 = 
2,e 1 = -l,c 2 = l) 


R = *J / 1 + iA^q™ 6 + iA 2 p^ s - A 1 A 2 M 12 (p 1 p 2 ) n (q 1 q 2 ) mS \ 

" 2 ° S \1 + iB lP ™q™ 5 + iB 2 p%q™ s - BiB 2 M 12 (p 1 p 2 )"(g 1 <fc)'"V ' 

However one can obtain a different nonlinear form by taking q n = f n /g n and 
r n = F n /G n , namely: 

__ Zif n+ i + Z 2 r n - 1 _ _ Ziq n+1 + Z 2 q n - i 

Qn — Qn * ~ , ( ^n — f'n * ~ , * 

Z 2 r „ + 1 + ^r n _i 2 2 ? ’n+l + 

where zi = 1 + + <5(1 + a 2 ), z 2 = 1 — icc — 5(1 + a 2 ). In the case of q n = r n 

the system becomes the Nijhoff-Capel lattice mKdV equation m but with 
complex coefficients. Of course a complete description imposes finding the Lax 
pairs, hamiltonian structure using r-matrices and so on, |18j . 

3 Second discretisation 

We may discretise the bilinear equations (4)-(7) in a simpler way. Namely we 
only discretise the dispersion equations (4) and (5) which depend explicitly on 


time and leave (6) and (7) unmodified. We obtain: 

fngn ~ fn'gn = 5(1 + « 2 ) {F n+1 G n -l ~ G n+1 F n _ x ) (20) 

F n G n — F n G n = 5(1 + a 2 )(f n +\g n -i — g n +\fn-i) (21) 

(1 + *a)i 7 ' IJ ,+iG ri _i + (1 — ia)F n -iG n +i = 2 f n g n (22) 

(1 + ia)f n +ig n -i + (1 - ia)f n -ig n +i = 2F n G n (23) 

This bilinear system is again completely integrable. It has TV-soliton solution of 
the same form, 

( N N \ 

( 24 ) 
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The dispersion relation and the phase factors are: 
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The nonlinear form can be easily recovered although the system now is more 
complicated and it will involve two auxiliary functions w n and v n . We divide 
(20) by g n g n and (21) by G n G n . Calling x n = fn/g n ,yn = F n /G n ,w n = 
Gn—iG n -\.\ /g n g n ,v n — f/n—i 9n-\-i / G n G rl we get: 

x n -x n = (5(1 + a 2 )(y n+1 - y n -i)w n 

y n -y n = <5(1 + a 2 )(x n+1 - x n -i)v n 

But one can see immediately that: 
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The factors in the parantheses can be computed easily from (22) and (23) by 
dividing them to G n -iG n +i and g n -ig n +i. Finally the nonlinear form of our 
sistem is: 

x n -x n = (5(1 + a 2 )(y n+1 - y n -i)w n 
Vn ~y n = < 5(1 + a 2 )(x n+1 - x n -i)v n 
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We can eliminate the auxiliary functions and z; n and we get the following 
higher order system: 


Xn+1 ^n+l _ Vn Vn X n -\-\X n -\-\ (1 d - ia) d - X n -\X n -\-\ (f ZGp 

Vn+ 2 -Un x n+ i-x n -i y n y n+2 {l + ia) + y n y n (l - ia) 

y n +i-y n +i x n -x n y n+ iy n+ i(l + ia) + y n -iy n+ i(l - ia) 


(28) 


(29) 


Zn+2 - X n Vn+1 - Vn-1 X n X n + 2 (l -f Z«) d" X n X„(l - ia) 

An important remark is that we have, the same phase factors and interaction 
term as in the differential-difference case, namely the system (4)-(7)(and that 
is why we kept the same notation). Also in both equations (28), (29) the step 
of time discretisation 5 dissapeared. So there is no trace of discretisation in the 
solutions, except the dispersion relation. This means that the structure of the 
soliton solution is the same at the level of tau functions. However, the nonlinear 
form is different. For instance the one and two-soliton solutions, 


x n = 


1 + PiPiQi 
1 + 7 i PiQ? 5 ’ 


Vn = 
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i + I’lPiq? 5 ’ 


mS 


%n — 
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are complex functions (since the phase factors are complex). But if we define 
new fields <j> n and if n by x n = exp(z<^„) and y n = exp(zz/> ra ) then the soliton 
solutions expressed by <f> n and i/j n are real functions and have the same shape 
as the ones in the first discretisation. 

Remark: This type of discretization involving auxiliary functions has been 
done also by Hirota and Tsujimoto in EMU, E3- In their examples also the 
form and phase factors of soliton solutions remain the same. We believe that 
this fact happens in more general cases, although must be proved rigurously. 

Also we could have discretised the system (1) directly by taking Q n (t) = 
T n (t)/$ n (t),R n (t) = H n {t)/T n (t) which give the bilinear system: 

DfTn • d’r, = H n +lT n —l — Hn-lT, 




DtH n • T n — r n _i^> n+ i 

+ C\T n <^ n + C2T^ = T n _iT n+ i 
Co Tn + ClH n T n + C 2 H^ = 
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Discretising the first two equations and after making the same steps as in the 
previous case we would have obtained almost a similar form: 


Qn Qn — ^(-^n+1 Rn— l)V n 
Rn Rn — ^(Qn+1 Qn— l)^kn 

ttt rr CO + ClQ n+ l + C 2 Ql +1 

™n+l *n r> \ r>2 

Co + C\R n + C 2 -/t^ 

CO + Cli?„+1 + C 2 ^J +1 

K+l = W n --- - - Ty2~ 

C 0 + CiQ n + c 2 Qi 

The main drawback of this system is that the soliton solution has a very com¬ 
plicated form. 


4 A new form of the Lotka-Volterra equation 


Let us give an nice and simple example related to the above construction, namely 
differential-difference Lotka-Volterra equation: 

du n , . 

— ^nf^n+l ^n—1 j 

Of course this equation can be seen as a particular case of the system (1) for 
Qn = R n = u n and Co = C 2 = 0, Ci = 1. However because the simplified 
form (2) has been obtained only for C 2 ^ 0 this equation cannot be studied as 
a particular case of the bilinear system (4)-(7). We consider the substitution 
u n = G n /F n and we get the following: 

D t G n ■ F n = ( G n +\F n -i — G n -iF n+ i) 


G F = F ,, F i 

y - T n ± n ± n+l x n— 1 

Now we are going to discretise only the first bilinear equation (imposing gauge- 
invariance in the right hand side) and the second one will remain the same. We 
shall obtain (t —> m6): 


G n F n — G n F n — 5(G n +\F n -i — Gn-iFn+i) 
G F = F , i F i 

n n+1- 1 - n— 1 


(30) 

(31) 


This bilinear system is integrable because it is equivalent (if we take G n = 
fn-ifn+ 2 , Fn = fnfn+i) with a quadrilinear equation reducible to: 


fnfn+l + Sf n -J n+2 - (1 + S)fJ n+1 = 0 (32) 


which is the integrable bilinear form of the discrete Lotka-Volterra (with u n = 
fn-lfn+ 2 /fnfn+l) 0 ] 

1 — 6 + Su n -i 


'1 - <5 + Su n+ i 
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However the nonlinear form recovered from (30) and (31) is different and we 
proceed as in the case of the system (28)-(29). Calling x n = G n /F n we get: 

Xn X n — ^(*£n+l X n —\)^Xn (33) 


w n+ i=io n ^±i (34) 

Xn 

This is a different variant of lattice Lotka-Volterra equation. We can solve the 
second linear equation in w n and we get: 


Xn X n — ${Xn +1 X n — l) 


k=—oo 


Xk 
Xk — 1 


or we can eliminate w n from the first equation and we find the following nice 
form: 

*£n+l X n -\-l X n -\-l *^77 — 1 X n ^ 

Xn X n X n -\-2 X n X n -\-\ 

The first two factors in the left hand side look like a discrete Schwartzian deriva¬ 
tive |T3] although the expression is not a cross ratio of four points but a cross 
ratio of diagonals of the two adjacent parallelograms formed by the six points. 
We think that this is a new equation although its symmetric form may be related 
in a way (unknown to us) to some well known one. 


5 Conclusions 

In this paper we have presented two integrable discretisations of a general bi¬ 
component differential-difference Volterra system. The main procedure was 
discretising differential Hirota bilinear operator and then recovery of nonlin¬ 
ear form with the aid of some auxiliary functions. This approach may lead 
to higher order nonlinear equations. We apply this procedure also to the well 
known Lotka-Volterra equation and we found a new discrete form However we 
started from a different bilinearization involving two tau functions. Relying on 
the fact that the structure of soliton solutions remains the same we believe that 
this procedure will be effective in discretising even nonintegrable equations of 
reaction-diffusion type (because it keeps the same structure of travelling waves). 
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Research 


6 Appendix: Proof of the N-soliton solution 

As we said in the introduction checking the existence of three soliton solution for 
a bilinear system is enough to prove integrability. However because the structure 
of the soliton solution is rather complicated we are going to sketch the proof 
for the N-soliton solution in the case of the first discretisation (for the second 


10 





one the phase factors and bilinear equations are simpler and everything goes in 
the same way). We follow the same procedure as in the old papers of Hirota 
he he m, he For simplicity lets call A = <5(1 + a 2 ) and redefine </,; —> re¬ 
plugging the N-soliton solution (12)-(15) into the first bilinear equation (8) we 
get 
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i<j 

where cimplies that the summation over families of indices p and // should 
be done under the requirements: 


Hi + n'i = M = 1, — >^0 
Hi = n'i = M = ^0 + 1, 

Hi = Hi = 0, i = v\ + 1, N 

Substituting the expressions of dispersion relation and phase factors (16)-(18) 
and introducing the multiindex a = fi — fj! we find that F = const.F where 


^0 


^0 


F = ^2 ( If (1 + ti^PV) ~ Ifi 1 + F^-Pi ai ) P (*«(-1 + cosh(cr i /c i )) - e i sinh(cr i fc l )(l + A)) - 


cr—±1 \i= 1 

/ ^0 


i= 1 


2=1 
V ^0 


-A ( JJ(1 + e*A Pi*)Pi* - JJ(1 + eiApi <Ti )p i ** J (zae»( 1 - c»cosh(a*£;;)) + sinh((Tifei)(l + A))} 

2=1 / 2=1 

x IT ( 6 * e j — COsh((T^ — (Tjkj )) 2 


^2=1 




Making the notation a:, = (e^) 1 / 2 , 2 = (1 + A + ia)/2, z* = (1 + A — ia)/2 the 
above relation becomes proportional to the following expression: 


p= em !!(-*«- 
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Z Xj + ZXj 3 )- 
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Now in order to show that (35) holds we have to prove that F = 0 for any 
n = 1, N. We shall prove this by induction. First, it is easily seen that this 
expression has the following properties: 

i) .F(:ri, x n ) is a symmetric and even function of x\,x n and invariant under 
the transformation Xj —> 1 / Xj for any j. 

ii) F{x i, x n )\ Xl -i = 0 

iii) F(xi, ...,x n )\ Xl=X2 = 2{l+/S.x\){l+/\x^ 2 ){-ia~z*x\+zx^ 2 ){-ia~z*x^ 2 + 
zx\)F{x 3 ,..., x n ). 

Now we assume that F = 0 for n— 1 and n— 2. From ii) and iii) F(xi,x n )\ Xi —i = 
0 and F(x i, x n )\ Xi - Xj = 0 for any i,j. Then, based on the property i), Hirota 
proved EE EE EZ3 that the expression of F can be factored by a function: 


nr=i(»? -1) 2 nS-(^ - ^) 2 (^ - x ) 5 


n?=i*?n 


(") r 4 r 4 
i<3 X i X 3 


which has 2 n + 2 n(n — 1) zeros of order 2. Accordingly F would have at least 
4n 2 polynomials in the numerator. On the other hand the function: 


n (n) 

F(xx n )Y[x 2 Y[xfx^ 

i=l i<j 


are polynomials of degree 2 n 2 + 6n at most. Hence F = 0 for any n. In the 
same way one can prove that the N-soliton solution obeys the other bilinear 
equations. 

For instance in the case of the third bilinear equation (10) introducing the 
expressing of N-soliton solution (12)-(15) we obtain 


E 

//,//= 0,1 


N 

n« 

3 =i 




N N / \ N ( N 

II«T + II-II a 7 b 7 II* 1 /n )3'; " j + 

U = 1 J =1 / i =1 v =1 


N 


J|(l -ia)q^ j pj j ^ 
f=i 


N 

x n n P^ i+ ^ n q^ i+ ^ )m = 0 (36) 

i<j 2=1 

From this expression we obtain by the same steps as above that the coefficient 
of the factor nT=] Pi ( li" YliLu +i Pi™Vi™ is proportional to the following expres¬ 
sion: 


k = E { 


(T —± 1 



— IOL ~ Z Xa 


— 2 <7j\ 

ZX 3 )- 
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^j(l + ia)( 1 + Ax 2 ^ )x 2 ° 3 + jj(l - ia)(l + Axj 2<T3 


u=i 


n< 

i<j 


j =1 


v^j -a.- 


^[(ia+2*x 2cr; '- 
j=i 


2<T ’ 2X7 2 " 3 )} 




Immediately one can see that K is an invariant, symmetric and even function 
satisfying the same properties K{x \,..., x n )\ Xl =i = 0, K{x \,..., x n )\ xl — X2 = 
2(1 + Axf )(1 + As( 2 )(-*a — z*x\ + zx^ 2 )(—ia — z*x (" 2 + zx\)K{x 3 , x„), so 
the induction method can be applied in the same way. 
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